Rapidity dependence of hadron production in central Au+Au collisions at y/s NN = 200 GeV 
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The rapidity and transverse momentum spectra for identified hadrons in central Au+Au collisions at -\fstm = 
200 GeV are computed in a quark combination model. The data of rapidity distributions for n ± , A"*, p(p) and 
net protons (p — ~p) are well described. We also predict rapidity distributions for K°, A( A), H~ (H + ) and Q.~ +Q + . 
The multiplicity ratios of charged antihadrons to hadrons as a function of rapidity are reproduced. The results 
for pr spectra of n ± , K*, p(p) and for the p/n ratios in a broader p T range agree well with the data. Finally the 
rapidity dependence of transverse momentum distributions for hadrons are given. 
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PACS numbers: 25.75.Dw, 25.75.-q 



I. INTRODUCTION 



The relativistic heavy ion collider (RHIC) at 
Brookhaven National Lab (BNL) provides a unique 
environment to search for the quark gluon plasma (QGP) 
predicted by lattice QCD calculations yj], and to study 
the properties of this matter at extremely high energy 
densities. A huge number of data have been accumulated 
and used to extract the information about the original 
partonic system and its space-time evolution. A variety 
of experimental facts from various aspects imply that 
the strongly coupled QGP has been probably produced 
in central Au+Au collisions at RHIC H d B 11 S 01 • 
Due to the confinement effects, one can only detect the 
hadrons freezed out from the partonic system rather than 
make direct detection of partons produced in collisions. 
Therefore, one of the important prerequisites for exploring 
the quark gluon plasma is the better understanding of 
hadronization mechanism in nucleus-nucleus collisions. 

The quark combination picture is successful in describ- 
ing many features of multi-particle production in high en- 
ergy collisions. The parton coalescence and recombination 
models have explained many highlights at RHIC, such as 
the high ratio of p/n at intermediate transverse momenta 
H, |sl [Toll and the quark number scaling of hadron elliptic 
flows [ 1 l,[l2|,[T3l]. Our quark combination model has been 
used to describe the charged particle pseudorapidity den- 
sities |3, hadron multiplicity ratios, pj spectra lfl5ll and 
elliptic flows II 1611 at midrapidity. 

Recently the STAR collaboration has measured the 
transverse spectra for identified baryons and mesons in a 
wider transverse momentum range in central Au+Au col- 
lisions at ys^7 = 200 GeV The BRAHMS collabo- 
ration has given the rapidity spectra Hill and the transverse 
momentum distributions for identified hadrons at different 
rapidities fl9l l20l |2lll . It provides a good opportunity to 
study the hadron production mechanism in a larger trans- 
verse momentum range and at different rapidities. In this 
paper, we use our quark combination model to study sys- 
tematically the rapidity and transverse momentum distri- 
butions for various hadrons in central Au+Au collisions at 
V^TaT = 200 GeV. 

The paper is organized as follows. In the next section we 



give a brief description of our quark combination model. 
In section III, the longitudinal and transverse distributions 
of partons in central Au+Au collisions at y/s^N = 200 
GeV are given. In section IV we compute rapidity dis- 
tributions of n ± , K 1 ", p(p) and the net proton multiplicity 
(p—p), and charged antiparticle-to -particle ratios as a func- 
tion of rapidity. We predict the rapidity distributions of K®, 
A(A), H~ (H ) and Q +Q + . In section V we present our re- 
sults for pr spectra of n ± , K* and p(p). We also give the 
p/n ratios in a larger pj range at y ~ and pj spectra of 
tt* and p(p) at other rapidities y ~ 1, 77 = 2.2 and y » 3.2. 
Section VI summaries our work. 



II. THE QUARK COMBINATION MODEL 

In this section we give a brief description of our quark 
combination model. The model was first proposed for high 
energy e + <?~ and pp collisions IL22L I23I I24I I25L L26L l27h . It 
has also been applied to the multi-parton systems in high 
energy e + e~ annihilations IT281 l29i l30l IJlll . Recently we 
have extended the model to ultra-relativistic heavy ion col- 
lisions suns. 

The average constituent quark number in nucleus- 
nucleus collisions can be written as (see e.g. Eq. (4) in 
Ref. H) 



{N q ) = 2[(a 2 + /? yW /2 - a]{N v ,J2) 



(1) 



where the two parameters are defined by (3 — l/2(V) and 
a = fim - j. Here (V) is the average inter-quark po- 
tential, m is the constituent quark mass given by m = 
(2m u + A s m s )/{2 + A s ), where m u and m s are the constituent 
masses of light and strange quarks respectively, and A s is 
the strangeness suppression factor. 

Our quark combination model describes the hadroniza- 
tion of initially produced ground state mesons (36 - plets) 
and baryons (56 - plets). In principle the model can 
also be applied to the production of excited states l24ll . 
These hadrons through combination of constituent quarks 
are then allowed to decay into the final state hadrons. We 
take into account the decay contributions of all resonances 
of 56 - plets baryons and 36 - plets mesons, and cover all 
available decay channels by using the decay program of 
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PYTHIA 6.1 (HI]. The main idea is to line up N q quarks 
and anti-quarks in a one-dimensional order in phase space, 
e.g. in rapidity, and let them combine into initial hadrons 
one by one following a combination rule. See section II of 
Ref. [15] for short description of such a rule. Of course, 
we also take into account the near correlation in transverse 
momentum by limiting the maximum transverse momen- 
tum difference Apr for quarks and antiquarks as they com- 
bine into hadrons. We note that it is very straightforward 
to define the combination in one dimensional phase space, 
but it is highly complicated to do it in two or three dimen- 
sional phase space [33]. The flavor SU(3) symmetry with 
strangeness suppression in the yields of initially produced 
hadrons is fulfilled in the model [ 22, 241. 
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HI. LONGITUDINAL AND TRANSVERSE 
DISTRIBUTIONS OF PARTONS 

The BRAHMS Collaboration has measured the ratios 
of antiparticle-to-particle and net-baryon rapidity density 
(dN( B _gJdy) as a function of rapidity in central Au+Au 
collisions at -\fsjm = 200 GeV IT341, l35ll . The results show 
that the ratios of K/K + , ~pjp decrease obviously with 
increasing rapidity and the dN, B _ B Jdy appears a valley 
shape in midrapidity region. This suggests that the rapidity 
distribution of net-quarks is different from that of newborn 
quarks. In this paper, we describe the momentum distribu- 
tions of newborn quarks by using a thermal phenomenol- 
ogy with a nonuniform collective flow. The rapidity distri- 
bution of net-quarks is given by three-sources relativistic 
diffusion model SIH. 

A. momentum distributions of newborn quarks 

We start with the momentum spectrum of quarks radi- 
ated by a stationary thermal source with temperature T 

E dW = d 3 n lh BcEe - E /T (2) 
J 3 p dy p-rdpjdip 

In the following text we give the spectra in terms of rapid- 
ity y = tanh^ 1 (pi/E) and transverse momentum px- 

The rapidity distribution of thermal quarks can be given 
by integrating Eq. (O over the transverse components l39ll 

dn t h T 3 lm 2 m 2 2 \ / m \ 

dy (27r) 2 w 2 T co shy cosh 2 y) \ T I 

(3) 

where m is taken to be constituent quark mass, i.e., m = 
0.34 GeV for light quarks/antiquarks and m = 0.5 GeV for 
strange quarks/antiquarks. Temperature of thermal source 
is taken to be T = 170 MeV, consistent with the phase 
transition temperature (T ~ 165 - 185 MeV ) from lattice 
QCD calculations. 

The momentum distribution of final-state hadrons mea- 
sured is anisotropic in high energy collisions. It rather has 



FIG. 1: (Color online) Rapidity spectra of quarks and antiquarks 
at hadronization in central Au+Au collisions at -\/s NN = 200 
GeV. 

TABLE I: Normalized expansion functions w(y) for light and 
strange quarks and antiquarks in central Au+Au collisions at 
= 200 GeV. 



u,umdd,d w(y) = 0.2exp( - |y| 2 - ffi / 14.85) 



s , s w(y) = 0.184 exp(-|>f 4 / 14.85) 



imprinted on the direction of the colliding nuclei. The 
boost-invariant longitudinal expansion scenario IT381 l39ll 
has explained such an anisotropy in terms of a boost- 
invariant longitudinal flow of matter with locally thermal- 
ized distribution. We apply the longitudinal flow to the 
parton level by adding up the spectra of individual ther- 
mal sources which are distributed in a rapidity region 
[y' min ,y' max ]- Similar to nonuniform longitudinal flow at 
hadron level 14011 . we speculate that the individual ther- 
mal sources are nonuniformly distributed in the rapidity 
region. The simplest way to describe this nonuniformity is 
to introduce a phenomenological expansion function w(y) 
in the following equation 

>™.v 

^(y) = N fdy'w(y')^(y-y'). (4) 

y',„i„ 

The transverse momentum spectrum of quark is not af- 
fected by this operation [39]. The integrity interval is cho- 
sen to be y' max = -y' min = y beam . N is the normalization 
constant. In the quark combination model, excluding the 
influence of resonance decays from final-state n + and K + 
rapidity distribution, we get the initially produced n + and 
K + rapidity spectra. We further inversely extract the ex- 
pansion functions w(y) for light and strange quarks and 
antiquarks, and obtain rapidity distributions of light and 
strange quarks and antiquarks. The results are shown in 
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Tableland Fig.[T]respectively. 

The transverse momentum spectrum of thermal quarks 
can be got by integrating Eq. (01 over rapidity using the 
modified Bessel function K\ 113911 



dn,. 



2n pT<ipT 



(5) 



dn t h 
2n pjt/px 



N 



which behaves asymptotically like a decreasing exponen- 
tial expi-mr/T), here transverse mass mj - yjp 2 + m 2 . 

Taken into account transverse flow of thermal source, 
the partons transversely boost by a flow velocity pro- 
file fi r {r) as a function of transverse radial positions r. 
/3,-(r) is parameterized by the surface velocity f3 s : /3,-{r) = 
[5 S (r/R max ) n . The transverse momentum spectrum of ex- 
pansion thermal source can be described as [39] 

r R "" I p T sink p\ jm T cosh p\ 

Jo >-dr mT I (-^)K l (-^—), 

(6) 

where Io is modified Bessel function, and p = tanhr l fi r is 
transverse rapidity. N is the normalization constant. The 
value of R max is taken to be R max - 13fm [41]. Values 
of parameters n, fl s for light and strange quarks and anti- 
quarks are extracted from the transverse momentum distri- 
butions of 7T° and K° s at midrapidity in our quark combina- 
tion model. 

Partons at high transverse momenta are mainly from the 
minijets created in initial hard collisions among nucleons. 
The transverse momentum distribution of minijet partons 
in the midrapidity region can be parameterized as [ 10, 42] 



dn 



2n Pt^Pt 



\B + p T ) 



The transition point (po) from thermal partons to mini- 
jet partons in p T spectra of quarks and antiquarks is deter- 
mined by the cross point of the two functions in Eq. (O and 
Eq. (O. Thermal partons dominate the transverse momen- 
tum below po while minijet partons dominate the trans- 
verse momentum greater than po. We extract the px spectra 
of quarks and antiquarks from measured tt° and K® spec- 
tra, then transition point po = 1.8 GeV/c is taken. For 
thermal quarks, We take n = 0.5 for light and strange 
quarks and antiquarks, B s = 0.45c, 0.56c for light and 
strange quarks/antiquarks respectively. Values of param- 
eters A, B and C for minijet partons are given in Ta- 
ble. |n] The scattering of minijet partons with thermal par- 
tons may lead to some change of the spectra of minijet par- 
tons and the smooth spectra around po. Normalized trans- 
verse momentum distributions for light and strange quarks 
and antiquarks at midrapidity in central Au+Au collisions 
at ysiviv = 200 GeV are in Fig. [2] 



B. rapidity distribution of net-quarks 

In this subsection, we use the three-sources relativistic 
diffusion model (RDM) JH, H3l to describe the rapidity 
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FIG. 2: (Color online) Transverse momentum spectra of quarks 
and antiquarks at midrapidity in central Au+Au collisions at 
-\[snn = 200 GeV. Minijet quarks have transverse momenta 
px > 1.8 GeV, while thermal quarks have transverse momenta 
p T < 1.8 GeV. 



TABLE II: Parameters for minijet parton distributions given in 
Eq.Q at midrapidity in central Au+Au collisions at yjs NN = 200 
GeV. 





A[l/GeV 2 ] 


B[l/GeV] 


C 


u, u and d, d 


2.32 


1.5 


7.25 


s , 1 


4.51 


1.5 


9.5 



(7) distribution of net-quarks. The time evolution of distribu- 
tion functions is given by a Fokker-Planck equation (FPE) 
in rapidity space |37j Hl+l |45l |4rjl 

j f [R(y, Of = ~[j(y)[R(y, r)f] + D y ^[R(y, t)f, (8) 



'8y 2 



where J(y) is drift function and D v is rapidity diffusion 
coefficient. 

Here, we use q = v = 1 corresponding to the standard 
FPE, and a linear drift function [37] 



J(y) = (y eq - y)/T y , 



(9) 



where r v is rapidity relaxation time, y eq is rapidity equilib- 
rium value and taken to be y eq = for symmetric systems. 

For the linear case above, exact solutions of standard 
FPE originating from R\^{y, t = 0) = 6(y + ybeam) can be 
obtained 04711 



l 



^o~\~(t) 



exp 



(y-<yi,2(0>r 
2a\ z {t) 



\ (10) 



where 



<yi,2(0> = +y beam exp (-t/T y ) 

o 2 12 (t)=D y T y [l-exp(-2t/Ty)l 



(11) 
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The available net-proton data at AGS ( Au+Au colli- 
sions at y*AW = 5 GeV) and SPS (Pb+Pb collisions at 
/saw = 17 GeV) are reproduced well with a superposition 
of two sources R\^,(y, in RDM. But it has been proposed 
in Ref. B43I1 that an expanding midrapidity source R 3 (y) 
emerges at RHIC. With this conjecture, the rapidity distri- 
bution of net-quarks in Au+Au collisions at ysAw = 200 
GeV can be written as 11361,1371 14311 

dN(y, t = Tim) 



dy 



= NiRi(y, t,„,) + N 2 R 2 (y, t,-„,) + N mid R 3 (y) 

(12) 



with the total number of net-quarks N\ + N 2 + N,„id equal 
to 3N part . Here t,„, is interaction time. The first and sec- 
ond items give the prediction of rapidity distribution for 
net-quarks at large rapidity, while the last R 3 (y) describes 
the rapidity distribution of net-quarks in the midrapidity 
region. 

There are only two adjustable parameters Ti nt /T y and 
D y Ty in R\,2(y,Tint)- Ti n t/Ty represents the diffusive time 
evolution of net-quarks; D y T y accounts for the widening of 
the rapidity distributions due to interactions among partons 
and partonic creations. In the present paper, D y T y has in- 
cluded the enhancement factor g( -\fs) caused by memory 
and collective effects |H3, BH H . In central Au+Au colli- 
sions at yfJffN = 200 GeV, values for parameters are taken 
to be t,„,/t v = 0.263 and D y T y = 0.9 respectively. 

The stationary solution of Eqj8]is approximate to Gaus- 
sian distribution ll47ll . It is regarded as R 3 (y) to describes 
the midrapidity valley for net-baryon rapidity distribution 
in Ref.Jla, H3, HH. Instead, for a better description of ra- 
pidity distribution of net-quarks, we obtain R 3 (y) by fitting 
the data of net-baryon in our model. The normalized poly- 
nomial fit for R 3 (y) in a rapidity region -y beam < y < y heam 
is 



R 3 (y) = 0.0715376 - 0.0189521/ - 0.004336/ 
+ 0.0000424407/ + 0.01001 M + 0.0194476|v| 3 . 



(13) 



The rapidity distribution of net-quarks in central Au+Au 
collisions at ysjviv = 200 GeV is shown in Fig. [3] Values 
of Ni,Nz,N m id are taken to be iVi = N2 = 214 and N mid 
607 for central Au+Au collisions at J% = 200 GeV. 



IV. RAPIDITY DISTRIBUTIONS OF IDENTIFIED 
HADRONS 

In this section, we use the quark combination model to 
compute the rapidity distributions of identified hadrons in 
central Au+Au collisions at -\fsm - 200 GeV. The to- 
tal number of constituent quark is given in Eq. CD. The 
values of parameters A s and /3 in Eq. CD ar e taken to be 
A s = 0.55, ft = 0.36 GeV" 1 respectivelylM [lj]. The ra- 
pidity distributions of newborn quarks and net-quarks are 
given respectively in Eq. © and Eq. (1121 1. The number of 
net-quarks is 3N par t and all the net-quarks are allowed to 




FIG. 3: (Color online) Rapidity distribution of net-quarks in 
central Au+Au collisions at yljviv = 200 GeV. The spectrum 
is multiplied by 1/3 to compare with the net-baryon data from 
BRAHMS Collaboration HI. 



take part in combination. Minijet partons are indispens- 
able to hadron productions at large transverse momenta. 
But from the transverse momentum spectra of partons at 
midrapidity in Fig. [2] we can see that the amount of mini- 
jet partons is about two order of magnitude smaller than 
that of thermal partons. The contribution from minijet par- 
tons to the yields and rapidity distributions of hadrons is 
negligible. Therefore, we don't distinguish between mini- 
jet and thermal partons in rapidity. With this input, we can 
calculate the rapidity distributions of various hadrons. 
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FIG. 4: (Color online) Pions and kaons rapidity densities as a 
function of rapidity in central Au+Au collisions at yljviv = 200 
GeV. The kaons yields were multiplied by 4 for clarity. The 
solid, dashed, dotted-dashed, dotted lines are our results for n + , 
n~ , K + and K respectively. The experimental data are given by 
BRAHMS Collaboration QJ]. 

In FigJU we show the rapidity distributions of charged 
pions and kaons in central Au+Au collisions at yjs NN = 
200 GeV. The pion yields are collected excluding the con- 
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tribution of hyperon (A) and neutral kaon decays. In 
quark combination picture, the existence of net-quarks will 
lead to the excess of direct produced K + over K yields, 
and decay contributions from other particles (K*°, <f>, Q) 
hardly make change to this excess. Therefore the excess of 
K + over K yields appears in the full rapidity range. While 
the directly produced n + and n~ yields can not reflect the 
influence of net quarks. Even if the decay contributions 
from other hadrons are considered, the difference between 
yields of n + and if~ is also very small compared with the 
large yields of pions. Thus the rapidity densities of n + and 
7r~ are nearly equal within the entire rapidity region. One 
can see that our model well describes the rapidity densities 
dN/dy of pions and kaons in the rapidity range covered for 
central Au+Au collisions at y/sNN = 200 GeV. 

The rapidity distribution of net-proton can reflect the 
energy loss of colliding nuclei and the degree of nuclear 
stopping. We calculate the rapidity distributions of proton, 
antiproton, and net-proton in central Au+Au collisions at 
-\[SffN = 200 GeV. The results are shown in Fig|5] No 
weak decay correction has been applied. The existence of 
net quarks also leads to the excess of directly produced 
proton over antiproton yields. Taken into account the de- 
cay contributions from hyperons, this excess will further 
increase. Our quark combination model can also success- 
fully explain the data of proton, antiproton, and net-proton 
in the rapidity region < y < 3. 
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FIG. 5: (Color online) Proton, antiproton and net-proton rapidity 
densities dN/dy as a function of rapidity in central Au+Au colli- 
sions at -\fsfm = 200 GeV. The solid and dashed lines in (a) are 
our results for proton, antiproton respectively, the solid line in (b) 
is our result for net-proton. The spectrum of net-quarks scaled is 
compared with that of net-proton. The errors shown with caps in- 
clude both statistical and systematic. The experimental data are 
given by BRAHMS Collaboration flUl. 

The process of quarks combination is a conversion 
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FIG. 6: (Color online) Antiparticle-to-particle ratios as a func- 
tion of rapidity in central Au+Au collisions at -yjstm — 200 
GeV. The dotted, dotted-dashed and solid lines are our results 
for rc fn* , K~ /K + , ~p/p respectively. The experimental data are 
from BRAHMS Collaboration lH. 



process of net-quarks to net-baryons simultaneously. In 
Fig(5fb), the computed rapidity distribution of net-proton 
nearly agrees with that of net-quarks except a small change 
of peak position around y » 4. This is mainly because 
of the mass effects in the processes of quarks combina- 
tion and resonances decays (see Ref. lfl4ll for explanation 
of this mass effects). Other baryons also have the simi- 
lar character. On the other hand, the rapidity distributions 
of hadrons are only slightly affected by rescattering in late 
stage of collisions [50]. Therefore, if quark combination 
mechanism is also applicable to large rapidity region, the 
measured net-baryon (B - B) rapidity distribution will re- 
flects quite truly the rapidity distribution of net-quarks just 
before hadronization. This provides a very beneficial con- 
dition to study the energy loss of colliding nuclei and the 
degree of nuclear stopping in nucleus-nucleus collisions at 
RHIC. 

Rapidity dependence of the ratios for antiparticles to 
particles are significant indicators of the dynamics of high 
energy nucleus-nucleus collisions l5Ul52ll . Since we have 
computed the rapidity spectra of charged pions, kaons, 
proton and antiproton, it is convenient to compute the ra- 
tios of yields of antiparticles to particles varied with rapid- 
ity. In Fig|6] we show the jr~/n + , K~/K + and ~p/p ratios 
as a function of rapidity in central Au+Au collisions at 
ys/VAi = 200 GeV. No weak decay correction has been ap- 
plied. The ratio of tt~/7t + is consistent with unity over the 
entire rapidity range, while the ratios of K~/K + and ~pjp 
decrease with increasing rapidity due to the influence of 
net quarks. Our results agree well with the data. 

From the above results, one can see that our model can 
well describe the rapidity distributions of primary charged 
hadrons. Among hadrons combined by light and strange 
quarks and antiquarks, the productions of hyperons S and 
O are less influenced by the decays of other hadrons, thus 
their rapidity distributions can reflect more directly the 
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FIG. 7: (Color online) Rapidity distributions for neutral kaon 
andhyperons A(A), E~ (H ), Q. +S1* in central Au+Au collisions 
at = 200 GeV. 

TABLE III: Rapidity densities dN/dy for neutral kaon K® and hy- 
perons A(A), H~ (E ), Q. +C1 at midrapidity in central Au+Au 
collisions at -y/s/wv = 200 GeV. The value of K° is the fit result 
according to its transverse momentum spectrum given by STAR 
Collaboration l53ll . other data are from STAR Collaboration 1 54] . 
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hadronization mechanism. We predict the rapidity spectra 
of neutral kaon and hyperons A(A), H~ (H ), Q +£2 + in 
central Au+Au collisions at y*MV = 200 GeV. The results 
are shown in Fig|7] We also compute the rapidity densities 
of these hadrons in midrapidity region and compare the re- 
sults with the experimental data. The data and calculation 
results are in Table Hill 



V. RAPIDITY DEPENDENCE OF TRANSVERSE 
MOMENTUM SPECTRA FOR IDENTIFIED HADRONS 

In this section, we calculate pj spectra of pions, kaons, 
protons and the p/n ratios in a larger pj range at y ~ in 
central Au+Au collisions at ~<Js NN = 200 GeV. We further 
compute the transverse momentum distributions of identi- 



fied hadrons at other rapidities y ~ 1, rj — 2.2 and y a 3.2. 

The transverse momentum distributions of newborn 
quarks and antiquarks have been given in section HI. Based 
on the experimental fact that n~/n + and ~p/p ratios at 
midrapidity are almost constants within a broad transverse 
momentum region {Vn\, the transverse momentum distri- 
bution of net-quarks is taken the same as that of newborn 
light quarks at midrapidity. We use Apr = 0.66 GeV for 
mesons and Apr = 1 .0 GeV for baryons respectively. With 
this input, we can compute the transverse momentum dis- 
tributions of various hadrons. 



A. Transverse momentum distributions of identified 
hadrons in a wider p T range at v ~ 

Identified baryon and meson distributions at large trans- 
verse momenta are available in Au+Au collisions at 
■\[snn = 200 GeV lfl7ll . It allows one to explore the particle 
production mechanism in a larger transverse momentum 
range. In this subsection, we compute the transverse spec- 
tra of identified hadrons at large p T and study the applica- 
ble extent of combination mechanism in a larger transverse 
momentum range. 

In FigJH] we show transverse momentum spectra of 
7t° and at midrapidity in central Au+Au collisions at 
s/snn — 200 GeV, which are used to determine the val- 
ues of parameters for parton px spectra. No decay correc- 
tions are applied to data. The transverse momentum spec- 
tra of pions, kaons, proton and antiproton are calculated at 
midrapidity in central Au+Au collisions at yfsiw = 200 
GeV. The results are in Figj9] The pions spectra are cor- 
rected to remove the feed-down contributions from K® and 
A. The feed-down contributions from A and E + are sub- 
tracted from the proton spectrum. Our quark combination 
model can also well explain the data in a large pj range. 

Recently the STAR Collaboration has measured the 
p/n + and ~pln~ ratios at large transverse momenta in cen- 
tral Au+Au collisions at yfstm - 200 GeV, and compared 
the data with the predictions of recombination and coa- 
lescence model 11711 . It is found that p/n + and p/n~ ratios 
peak at pj ~ 2-3 GeV with values close to unity, decrease 
with increasing px, and approach p/n + as 0.4, ~p/n~ w 0.25 
at px ^7 GeV. These models can qualitatively describe 
the pCp)/7T ratio at intermediate px but in general under- 
predict the results at high p T II 1711 . We also compute the 
pln + and p/n~ ratios in the range px < 12 GeV. The cal- 
culation results are shown in Fig[l0] Our results agree with 
the experimental data in the full p T range. 

One can see from the above results that the quark com- 
bination model can reproduce the productions of identified 
hadrons quite well in a larger transverse momentum range. 
It suggests that quark combination mechanism may still 
play an important role at high transverse momentum up to 
12 GeV/c. 
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FIG. 8: (Color online) Transverse momentum spectra of n° (a), 
(b) at midrapidity in central Au+Au collisions at sfstw = 200 
GeV. The solid lines are our results. The data of n° are given 
by PHENIX Collaboration 0], and data of K° are from STAR 
Collaboration |53fl. 



B. Transverse momentum spectra of hadrons at other 
rapidities y ~ 1, 77 = 2.2 and y ~ 3.2 

A variety of experimental facts at RHIC indicate that the 
QGP has been probably produced at midrapidity. Based 
on this indication, quark combination picture has suc- 
cessfully explained the hadron production at midrapidity 
lfl5i Ull However, can the QGP extend to forward 
rapidity? If it does, how much rapidity can it extend to? 
Can combination picture describe the hadron production 
at forward rapidity yet? Fortunately, the BRAHMS Col- 
laboration recently has measured the pj spectra of identi- 
fied hadrons at different rapidities [ 19, 20, 21]. It provides 
an opportunity to study the hadrons production at differ- 
ent rapidities and rapidity dependence of hadrons trans- 
verse momentum spectra. We have described the longi- 
tudinal distributions of identified hadrons in section IV. 
The transverse momentum spectra of hadrons at midrapid- 
ity are well reproduced in the above subsection. We now 
make a straightforward extension to the forward rapidity 
region, and this extension should be made with no change 



in the quarks px spectra before we draw some useful phys- 
ical information from it. 

We compute the transverse momentum spectra of pi- 
ons and protons at y ~ 1 in central Au+Au collisions 
at ysMV = 200 GeV, the results are shown in Fig fTTT a). 
The proton and antiproton spectra are corrected to remove 
the feed-down contributions from A and A weak decays. 
No significant changes for the transverse spectra of pi- 
ons and proton(antiproton) are observed within one unit 
around midrapidity. It suggests that hadrons production at 
this rapidity are still dominated by combination in inter- 
mediate transverse momentum range. The QGP still exist 
at y ~ 1 and it may extend much more than this rapidity. 

In Fig fTTT b). we show the transverse momentum spectra 
of pion and antiproton at 77 = 2.2(2.14 < 77 < 2.26) in 
central Au+Au collisions at yfsim = 200 GeV. The cor- 
rection for feed-down from the (anti-)lambda decays has 
been applied, whereas the contamination of pions due to 
weak decays has not been corrected. Our model can also 
well describe the hadronic production at this rapidity. 

We further compute the transverse spectra of pions and 
(anti-)proton at y a 3.2, the results are shown in Fig fTTT c). 
No correction for decay or feed-down has been applied to 
data. One can see that computation results are in good 
agreement with the data for pions, but can't reproduce the 
data of proton and antiproton. This is probably because 
the transverse momentum distribution of net-quarks varies 
with rapidity in nucleus-nucleus collisions at RHIC. The 
amount of net-quarks is marginal compared with that of 
newborn quarks in midrapidity region. Net-quarks may 
be highly thermalized due to interactions with newborn 
quarks. Its transverse momentum distribution will asymp- 
totically tend to that of newborn light quarks. However the 
amount of net-quarks can approach and even exceed that 
of newborn quarks at large rapidity. Therefore the ther- 
malization extent of net-quarks at large rapidity is much 
smaller than that at midrapidity. Its transverse momen- 
tum distribution at large rapidity may be obviously differ- 
ent from that of newborn quarks. Just as analyzed in sec- 
tion IV, the yields of proton and antiproton are more sensi- 
tive to the existence of net-quarks than that of pions. Thus 
the rapidity dependence of transverse momentum distribu- 
tion for net-quarks should be mainly embodied by that of 
proton and antiproton, rather than pions. It is one of the 
probable reason why our quark combination model can de- 
scribe well the data of pions but for proton and antiproton 
aty ~ 3.2. 

From the above results, we can see that the quark com- 
bination model is able to describe the transverse momen- 
tum distribution of identified hadrons in a broad rapidity 
range < y < 3.2. This implies that quark combination 
hadronization mechanism is still applicable to large rapid- 
ity (at least y - 3.2). On the other hand, the BRAHMS 
Collaboration has measured the nuclear modification fac- 
tor Raa for identified hadrons at forward rapidity in central 
Au+Au collisions at y/s NN = 200 GeV 112 ill , and the data 
also indicate the existence of quark combination mecha- 
nism at forward rapidity. The continued suppression seen 
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FIG. 9: (Color online) Transverse momentum spectra of pions (a), kaons (b), proton and antiproton (c) at midrapidity in central Au+Au 
collisions at -\/s NN = 200 GeV. The solid lines are our results for n + , K + and p, and the dotted-dashed lines are for n~, K~ and ~p 
respectively. The data of pions, proton and antiproton are from STAR Collaboration fT*7tl . and kaons are from BRAHMS Collaboration 
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FIG. 10: (Color online) The ratios of p/n + and ~p/n~ as a func- 
tion of transverse momentum at midrapidity in central Au+Au 
collisions at sjs NN = 200 GeV. The shaded boxes represent the 
systematic uncertainties. The solid lines are our results. Experi- 
mental data are given by STAR Collaboration fl7ll . 



ferent physics from that at midrapidity. The absence of 
suppression in particle production at large transverse mo- 
mentum in d+Au collisions shows that the suppression in 
central Au+Au collisions at midrapidity is not an initial- 
state effect but a final-state effect of the produced density 
medium (jet quenching) 11571. l58ll . However, the initial-state 
parton saturation effects are more evident at large rapid- 
ity i59l l60ll . which has been testified by the high rapid- 
ity suppression measured in d+Au collisions 16 ill . But the 
decreasing of the parton medium density at large rapidity 
will lead to the smaller jet energy loss, thus the less sup- 
pression caused by jet quenching can be expected. There- 
fore the suppression measured at large rapidity in central 
Au+Au collisions may be the result of a compromise be- 
tween initial- and final-state effects i62ll . Maybe, the same 
parton px spectra at y = 3.2 with that at y = used in our 
work is just a synthetic embodiment of the two effects. 



in the Raa at y = 3.2 further suggests that there is dif- 
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FIG. 1 1 : (Color online) Transverse momentum spectra of pions and protons at y ~ 1 (a), r\ = 2.2 (b) and y = 3.2 (c) in central Au+Au 
collisions at ^]s NN = 200 GeV. The solid lines are our results for 7r + and p, and the dotted-dashed lines are for 7r~and ~p. The data of 
panel (a),(b),(c) are given by BRAHMS Collaboration CUMIH- 



VI. SUMMARY 

Using the quark combination model, we study the rapid- 
ity and transverse momentum distributions for identified 
hadrons in central Au+Au collisions at -\fsm = 200 GeV, 
incorporating the collective expansion of the hot medium. 
We use the measured rapidity spectra of n + and K + to 
extract the values of parameters and the explicit form of 
the nonuniform expansion function. The rapidity distri- 
butions of net-quarks are described by a relativistic diffu- 
sion model i37ll . We compute the rapidity distributions of 
identified hadrons. The results agree with the experimen- 
tal data. The existence of net quarks leads to the excess 
of K + over K and that of proton over antiproton in full 
rapidity range. The difference between yields of n + and 
n~ resulting from net-quarks is very small compared to the 
large yields of pions. We further predict the rapidity spec- 
tra of K Q S , A(A), "Er (H + ) and + Q + in central Au+Au 
collisions at yfstm = 200 GeV. We also calculate the trans- 
verse spectra of identified hadrons and the ratios of pro- 
ton/ antiproton to pions in a large transverse momentum 



range at midrapidity. Our model reproduces the data quite 
well, which suggests that combination scenario may still 
play an important role in hadronization at high pj up to 12 
GeV/c. Finally we give rapidity dependence of transverse 
momentum spectra of hadrons. Using the same transverse 
momentum spectra for quarks, we calculate the transverse 
momentum spectra at various rapidities y ~ 0, 1, r\ = 2.2 
and y » 3.2. The results are in good agreement with the 
data except those for proton and antiproton at y « 3.2. 
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